Origin of ferroelectricity in high- Tc magnetic ferroelectric CuO 
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Cupric oxide is a unique magnetic ferroelectric material with a transition temperature significantly 
higher than the boiling point of liquid nitrogen. However, the mechanism of high-Tc multiferroicity 
in CuO remains puzzling. In this paper, we clarify the mechanism of high-Tc multiferroicity in CuO, 
using combined first-principles calculations and an effective Hamiltonian model. We find that CuO 
contains two magnetic sublattices, with strong intrasublattice interactions and weakly frustrated 
intersublattice interactions, which may represent one of the main reasons for the high ordering 
temperature of the compound. The weak spin frustration leads to incommensurate spin excitations 
that dramatically enhance the entropy of the mutliferroic phase and eventually stabilize that phase 
in CuO. 

PACS numbers: 75.85.-l-t, 71.20.-b, 75.25.-j 
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"Magnetic ferroelectric" materials, in which ferroelec- 
tricity is induced by magnetic ordering, have attracted 
intense interest [l|, Q ■ The strong magnetoelectric (ME) 
coupling in these materials opens a new path to the de- 
sign of multifunctional devices that allow the control of 
charges by the application of magnetic fields or spins 
through applied voltages. However, nearly all current 
magnetic ferroelectric materials are strongly frustrated 
magnets with very low ordering temperatures (~ 30 
- 40 K), several times lower than the temperatures ex- 
pected from the strengths of their spin interactions. Low 
critical temperature is one of the major factors that limit 
the application of these important materials. Therefore, 
a new mechanism that allows high-temperature magnetic 
ferroelectricity is highly desirable. 

Recently, CuO has been found to be multiferroic at 
Tc=2iQ K, which is much higher than the critical tem- 
peratures of any other magnetic ferroelectric materials 
[3|. CuO undergoes two successive magnetic phase tran- 
sitions when it is cooled from room temperature to a tem- 
perature near zero. Neutron scattering experiments Q 
have shown that at temperatures below Tjvi=213 K, the 
spin structure is coUinear antiferromagnetic (AFMl) [see 
Fig. [l^a)]. Between Tmi and TAr2=230 K, the spin struc- 
ture becomes non-collinear and slightly incommensurate 
(AFM2) [see Fig. IHb)], with a modulation vector of 
Q= (0.006, 0, 0.017). An electric polarization of 160 ^iC 
mT^ along the b axis also develops in the AFM2 phase. 
Elucidation of this unusual high-temperature multiferroic 
behavior may provide useful information in the search 
for novel room-temperature magnetic ferroelectric mate- 
rials. However, the mechanism that stabilizes the multi- 
ferroicity in CuO is still not understood, and it remains 
the subject of significant debate Q- The phase dia- 
gram of CuO has been studied via a phenomenological 
approach by Toledano et al. 0, which based on sym- 
metry considerations only, could not reveal the micro- 
scopic mechanism that stabilizes the multiferroic phase 




FIG. 1: Schematic sketch of the magnetic structures of (a) 
the colhnear AFMl phase, and (b) the noncolhnear AFM2 
phase. The black arrows, yellow crosses, and blue circles de- 
note the spin directions associated with Cu ions. The black 
circle in (a) indicates an inversion center, (c) A sketch of su- 
perexchange interactions Ji to J-j. The single lines, double 
lines, and dashed lines represent the three types of exchange 
interactions between Cu ions. 



in CuO. Giovannetti et al. performed Monte Carlo sim- 
ulations on an effective Hamiltonian model of CuO 
However, in the simulation the spins were artificially con- 
strained into only four possible directions, which dramat- 
ically changes the free energy of the system and misses 
important physics in CuO. 

In this work, we clarify the mechanism of high-Tc mul- 
tiferroicity in CuO using combined first-principles calcu- 
lations and simulations based on an effective Hamiltonian 
model. We find that CuO contains two magnetic sublat- 
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tices, with strong intrasublattice interactions and weakly 
frustrated intersublattice interactions, which might rep- 
resent one of the main reasons that the compound ex- 
hibits a high ordering temperature. The weak spin frus- 
tration leads to incommensurate spin excitations that 
dramatically enhance the entropy of the AFM2 phase and 
eventually stabilize the mutliferroic phase. This mecha- 
nism is novel and differs from previously proposed mech- 
anisms [1, @] . This work suggests that high- Tc magnetic 
ferroelectric materials can be sought in weakly frustrated 
magnets similar to CuO. 

The crystal structure of CuO is monoclinic and con- 
tains four chemical units per unit cell. The AFMl spin 
structure is composed of two antiferromagnetic (AFM) 
spin sublattices, shown in two different colors in Fig. [T] 
(a), in which Cu ions have the same h values in each sub- 
lattice. The spin chains along the [101] direction are anti- 
ferromagnetic, and are labeled chain I and chain II for the 
two sublattices, whereas the chains along the [101] direc- 
tion are ferromagnetic. In the AFMl phase, all spins are 
aligned in the h direction, whereas in the AFM2 phase, 
chain II rotates perpendicularly to chain I. 

We perform ab initio calculations on CuO, with 
non-collinear spin-polarized local density approximation 
(LSDA) implemented in the Vienna ab initio simulations 
package (VASP) 0,0. The on-site Coulomb interactions 
U~7.5 eV are included for Cu ions in a rotationally in- 
variant scheme Q . The spin-orbit coupling is considered 
in the calculation unless otherwise noticed. To accommo- 
date the spin structures, we use a 2x1x2 CuO supercell 
that contains 32 atoms. For the AFM2 structure, we ne- 
glect the small incommensurate component of the spin 
structure [i.e., we set Q= (0, 0, 0)], and rotate the spin 
directions of chain II 90°, so that it lies in the ac plane. 
The spins form cycloidal spirals along both the a and c 
axes in the AFM2 phase. The incommensurate compo- 
nent of the magnetic modulation vector Q is extremely 
small, and it should not affect the calculated electric po- 
larization because P (x Si x UM. More details of the 



first-principles calculations can be found in Ref.|ll|. 

We first determine the crystal structure of CuO un- 
der the AFMl spin configuration. The room tempera- 
ture crystal structure is monoclinic and of space group 
C2/c, with inversion symmetry. However, the magnetic 
structure of the AFMl phase only has P2i/c symmetry. 
Therefore, after relaxation, the crystal structure is also 
reduced to P2i /c symmetry because of "exchange stric- 
tion" effects [I4]. Further analyses {TIJ show that the 
Cu and O ions deviate from their high symmetry sites by 
approximately 10"'^ A. This distorted structure preserves 
the inversion symmetry; therefore it has no net polariza- 
tion. The inversion center is shown in Fig. [Ija). An 
inversion operation about the inversion center changes 
spin chain I to chain II. 

We then fix the spin orientations to the AFM2 configu- 
ration and relax the crystal structure again to obtain the 



crystal structure of the AFM2 phase. It turns out that 
the AFM2 configuration is a local minima of the total 
energy. The results change little if the spins are let free 
to rotate after the initial spin configuration is set to the 
AFM2 phase. The calculated total energy of the AFM2 
phase is greater than that of the AFMl phase by approx- 
imately 0.33 meV per atom, which is consistent with ex- 
periments that have shown the AFM2 phase to appear at 
a higher temperature than that at which the AFMl phase 
appears [3. \WL . Before structural relaxation, the total 
energy of the AFM2 phase is approximately 0.08 meV per 
atom higher than that of the AFMl phase. This energy 
difference is primarily attributed to the spin anisotropy 
energy, which can be observed when the spin-orbit cou- 
pling is turned off, which reduced the energy difference 
between the two phases to 0.03 meV per atom. The re- 
maining difference of 0.25 meV per atom exists because 
the ionic distortion in the AFMl phase is substantially 
larger than that in the AFM2 phase, as shown below. 

The AFM2 spin structure exhibits P2i space group 
symmetry, to which the rotation of spin chain II does 
not give inversion symmetry, as can be seen in Fig. [Ub). 
The crystal structure is distorted by the Dzyaloshinskii- 
Moriya (DM) interaction l^- l3l , which breaks the inver- 
sion symmetry. When compared with the high symmetry 
structure, all the oxygen ions are shifted in the +b direc- 
tion by approximately 7x10^^ A, whereas all Cu ions 
are shifted in the -b direction by a similar amount. The 
ionic distortion in the AFM2 phase is approximately two 
orders of magnitude smaller than that driven by the "ex- 
change striction" effects in the AFMl phase. 

Next, we calculate the electric polarization using the 
Berry-phase theory of polarization The calculated 
total polarization is approximately 90 /iC m~^ in the -b 
direction, which is somewhat smaller in magnitude than 
the experimental value of 160 /iC m^^ along the b axis. 
The agreement between the theoretical calculations and 
experimental values is reasonable, given that the current 
functionals are not adequate to treat the subtle correla- 
tion effects in magnetic ferroelectric materials . 

In magnetic ferroelctric materials, the transition tem- 
peratures are predominantly determined by the mag- 
netic exchange interactions [3|. We extract the superex- 
change interactions Js of CuO using a Heisenberg model, 
Hm — — J2ij Jij^i • — ■ ^i)^ from the calculated 

total energies of the different spin configurations in the 
symmetrized C2/c crystal structure with spin-orbit cou- 
pling. K is the anisotropic energy due to spin-orbit cou- 
pling. There are seven J values in total, which are shown 
in Fig. 1(c). Among them, Ji is the exchange interac- 
tion between nearest-neighbor Cu atoms along the [101] 
direction, J2 is the interaction between nearest-neighbor 
Cu atoms along the [101] direction, Jj is the interaction 
between the nearest-neighbor spins of the same sublat- 
tice along the b direction, and J3 and J4 are the inter- 
sublattice exchange interactions. The fitted value of Ji = 
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FIG. 2; The polar lattice distortions 112 as functions of tem- 
perature. The maximum value of the lattice distortion is set 
to 1. Blue curve: the result of the 12x12x12 lattice. Red 
curve: the result of the 48x4x48 lattice without frustrated 
spin interactions. Black curve: the result of the 12x12x12 
lattice, with the spins are constrained in 4 perpendicular di- 
rections. The lattice distortion values are amplified by factors 
of 20, or 2 as indicated in the figure. 



-51 meV is in good agreement with the value determined 
from neutron scattering experiments, J=67±20 meV[^. 
The fitted value of J2—8.6 meV, is only approximately 
1/6 the value of which is consistent with the quasi- 
ID model and the fitted value of J7=9.87 meV. 

The fitted inter-sublattice coupling value J3=4.9 meV 
and J4=7 meV are weak, because the Cu-O-Cu bond an- 
gles are close to 90° for these two Js [l^. In the AFMl 
and AFM2 phases, symmetry causes additional mutual 
cancellation of J3 and J4; therefore, the energy cost of 
chain II rotation is low, as was discussed in the previous 
paragraph. We also calculate the next-nearest-neighbor 
interactions, J5 and Je. We find that J5 and Je are highly 
asymmetric. Interaction J5 is significant, with a value of 
-12 meV, whereas Je is only 2.1 meV. These values agree 
well with those of Filippetti and Fiorentini 21 1. 

The intrasublattice interactions Ji, J2 and J7 essen- 
tially determine the ground state spin structure, AFMl. 
The J5 next-nearest-neighbor interactions further favors 
the antiferromagnetic spin chain along the [101] direction, 
whereas, Jq only adds a small frustration to this configu- 
ration. The major competing interactions are those of the 
inter-sublattice interactions, J3 and J4. The weak incom- 
mensurateness of the spin spiral caused by the frustrated 
exchange interactions J3 , J4 is consistent with that effect 
in this material, because the spin competition is small. 
We calculate the ordering temperature of the Heisenberg 
model by a Monte Carlo simulation, in which all exchange 
interactions are forced to be ferromagnetic and we ob- 



tain Tc= 311 K [2^. This temperature is only about 1.5 



times greater than the Tc value of the AFMl phase, (in 
RMn2 05 [iH. the ratio is approximately 3 - 4.), which 
also indicates that the spin frustration is weak in CuO. 
The lack of strong competing interactions in this com- 
pound may explain the high spin-ordering temperature 
of CuO. 



To study the mutliferroic phases of CuO, we simulate 
the full Hamiltonian model with spin-lattice interactions, 
H = i/ph + Hm + Hi + Hum, where, 

Hph = -E^o + X! ^miujlul{k) + ]-m2Ujlul{k) , (1) 



Hi = ^^4ui(fc)S,.S,, 



i.j k 



Hum = X! D[u2{i)h x ei^,+i] ■ (S., x S 
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Here Eq is the energy of the high-symmetry structure 
without magnetic interactions. Since the symmetry- 
lowering displacement is extremely small, we treat this 
displacement in the AFMl (AFM2)_phase as the nonpo- 
lar(polar) phonon mode ui {U2) [23[. Ui(fc) and U2(k) 
are the phonon modes in the kth unit cell, and mi (7712) 
and oji (012) are the reduced mass and frequency of the 
non-polar modes (polar modes), respectively. For sim- 
plicity, we neglect the phonon dispersion in the simu- 
lation, which has little effect on the results [2^. Hi is 
the isotropic spin-lattice interaction caused by exchange 
striction effects, J[j = ^jf- In CuO, only J3 and J4 are 
involved in the nonpolar lattice distortion, and all other 
J's are canceled out by symmetry. Hum is the DM inter- 
action term, which sums over the nearest neighbor spin I 
along the a and c axes, where b is the unit vector along 
the b (polar) axis and ei^^i is a unit vector connecting 
the i-th and the i + l-th Cu atoms. Unlike the model 
in Ref. @, our model explicitly includes the lattice de- 
gree of freedom, especially the non-polar modes, which 
are missing in the model of Ref. [5j . 

For simplicity, we redefine Ui (^2) to be a dimension- 
less parameters that takes the value of unity at the low- 
symmetry state of the AFMl (AFM2) phase, and we as- 
sign spin moments |Si|=1.0. All parameters in the model 
can be obtained by fitting to the total energies of the first- 
principles calculations [23|. We use J3=J4=0.8163 meV, 
which is approximately 1.5 times greater the values ob- 
tained by first-principles calculations . To obtain the 
correct ground state, we reduce the values to J3=2.45 
meV, and J4=3.5 meV, which are approximately half of 
the fitted values, miwf =0.4, m2U}j=0A9D. i:)c=0.8723 
meV is the critical value at which the energies of AFMl 
and AFM2 phases are degenerate. Our main results are 
valid in a reasonable range of the values of parameters 
for the model Hamiltonian, provided J3 and J4 are larger 
than O.27J3 (to ensure that the ground state is in the 
AFMl phase) and J3 and J4 are larger than half of the 
fitted values. 

We simulate the effective Hamiltonian model in the 
temperature range from 29 to 387 K, using a replica- 
exchange Monte Carlo method 0|. We perform the 
Monte Carlo simulations on the nxmxn lattices (71= 12 
- 72, m=2 - 12) with periodic boundary conditions. The 
electric polarization P = (M2). Typical simulation re- 
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FIG. 3: (a) The lattice distortions ui (black) and U2 (red) 
as functions of temperature. The maximum value of the lat- 
tice distortion is set to 1. (b) The transition temperature T2 
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FIG. 4: (a) The entropies as functions of the temperature of 
the 48x4x48 lattice with (black line) and without (red line) 
frustrated spin interactions, (b) The entropy difference be- 
tween the systems with/without frustrated spin interactions. 



suits for the n=12 lattice are shown in Fig. [51 There is 
no stable AFM2 phase, even with a large DM interaction, 
D=0.5Dc- If we artificially constrain the spin to only four 
possible directions at 90°, as was done by Giovannetti at 
al. 0, we obtain a weak AFM2 phase at approximately 
300 K, as shown in Fig. [5J that is similar to the one 
obtained by Giovannetti et al. The above results 

suggest that the DM interaction alone cannot stabilize 
the AFM2 phase in CuO. 

Surprisingly, if the simulation lattice size is increased 
to n >32, a well defined AFM2 phase is obtained. The 
simulation results of a 48x4x48 lattice with a small 
£'=0.066 Dc are presented in Fig. [SJa). A paraelec- 
tric (PE) to ferroelectric transition clearly occurs near 
r2=180 K. At Ti= 125 K, the polarization suddenly 
drops to zero, accompanied by the appearance of the 
non-polar lattice distortion, which indicates that the sys- 
tem enters the AFMl phase. This result is in excellent 
agreement with experiments. Fourier analyses of the spin 
structures also confirm the above results. A finite size 
scaling for T2 with m gives £2=192 K as 1/m approaches 
zero [see Fig. EJb)], which slightly underestimates the 
experimental value of T2=230 K i3]. 

It is initially surprising that the AFM2 phase is not sta- 



ble in the small simulation lattice but survives in the large 
simulation lattice. Fourier analyses of the spin structures 
suggest that, in the n=48 lattices, many incommensurate 
spin components are present, which are forbidden in the 
small (e.g., n=12) lattices. The incommensurate spin 
components are generally known to arise from the frus- 
trated spin interactions, which suggests that the small 
spin frustration might play an important role in stabiliz- 
ing the AFM2 phase. We then perform the simulations 
using the 48x4x48 lattice, but set the frustrated spin 
interactions to zero, i.e., J3=J4=Je=0. The results are 
shown in Fig. [2l Remarkably, the results are very simi- 
lar to those obtained for the n=12 lattice, and no stable 
AFM2 phase is found. We also carry out simulations 
by using the effective Hamiltonian and parameters given 
in Ref.[5| and obtain similar results. This result further 
confirms the assumption that the weak spin frustration 
is essential to stabilize the AFM2 phase, and contradicts 
to previous conclusions [1, @] . 

To understand why the weak incommensurate spin 
components can stabilize the AFM2 phase, we calculate 
the entropy of the system as a function of the temper- 
ature using a multihistogram reweighting technique [25|. 
The entropies of the 48x4x48 lattice with or without 
spin frustration are compared in Fig. |4l The entropy of 
the non-frustrated system increases smoothly as temper- 
ature increases from 120 K to 130 K. The entropy of the 
71=12 lattices with frustrated interactions (not shown) 
exhibits a similar temperature dependence. In contrast, 
an obvious entropy jump occurs near Ti=125 K in the 
48x4x48 lattice with spin frustration. This increase can 
be seen more clearly in Fig. |4] (b), which depicts the 
difference of the system entropies with and without spin 
frustrations in the temperature range 50 -160 K. There- 
fore, we conclude that the incommensurate spin excita- 
tions caused by the spin frustration greatly enhances the 
entropy of the system in the AFM2 phase and stabilizes 
the phase. 

To explore the role of the DM interaction in the phase 
transitions, we tune the DM interaction strength D from 
zero to Dc- In the range of I? < Dc, transitions from the 
paramagnetic-PE phase to the AFM2-ferroelectric phase 
to the AFMl-PE phase are always present. As D in- 
creases, the temperature range of the AFM2 phase also 
increases. When D > Dc, only the transition from the 
paramagnetic-PE phase to AFM2-ferroelectric phase ex- 
ist, and there is no AFMl phase Interestingly, when 
we set D—0, we still obtain a stable AFM2-like incom- 
mensurate phase, but it has no net electric polarization 
because there is no unique polarization axis in this case. 
The DM interaction that breaks the rotational symme- 
try of the nearly degenerate spin structures in the AFM2 
phase and generates the polarization axis |5|- However, 
in the case of a strong magnetic anisotropy, electric po- 
larization is also possible. 

Understanding the mechanism of the multiferroicity 
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in CuO provides important guidance in the search for 
new high-Tc magnetic ferroelectric materials. We found 
that CuO contains two magnetic sublattices, with strong 
intra-sublattice interactions and weakly frustrated inter- 
sublattice interactions, which may represent one of the 
main reasons that the compound has a high ordering tem- 
perature. Monte Carlo simulations suggest that the in- 
commensurate spin excitations, caused by the weak frus- 
trated interactions dramatically enhance the entropy of 
the multiferroic AFM2 phase, and they eventually stabi- 
lizes the phase. The DM interaction break the magnetic 
rotational symmetry, which causes lattice distortion and 
lead to electric polarization in the AFM2 phase. This 
mechanism is distinct from previously proposed mecha- 
nisms Qfor CuO. One of the fascinating features of 
magnetic ferroelectric materials is their rich phase dia- 
grams, which arises from the competing interactions in 
these materials. However, few reports on the mechanism 
of these phase transitions at the microscopic level have 
been published. The methods developed in this work may 
be useful in elucidating the complex magnetic ferroelec- 
tric phase transitions in general magnetic ferroelectric 
materials. 
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